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Cosmology with a three-form field interacting with cold dark matter is considered. In comparison 
to coupled scalar field quintessence, the new features include an effective pressure contribution to the 
■ field equations that manifests both in the background and perturbation level. The dynamics of the 

background is analyzed and new scaling solutions are found. A simple example model leading to a de 
Sitter expansion without a potential is studied. The Newtonian limit of cosmological perturbations 
£NJ ■ is derived and it is deduced that the coupling can be very tightly constrained by the large-scale 

^ ' structure data. This is demonstrated with numerical solutions for a model with nontrivial coupling 

, and a quadratic potential. 
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O ' I- INTRODUCTION 

u' 

The tiny magnitude of the cosmological constant in the standard ACDM model is not theoretically understood. In 
attempts to tackle this problem the cosmological constant A is often replaced by a dynamically evolving entity called 
dark energy as the agent behind the Universe's present acceleration The theory underpinning the inflationary 
paradigm that extends the standard big bang model to the earlier universe also remains to be established. Most often 
the field responsible for dark energy 0, El as well as inflation is a scalar field, but it is of interest to study possible 
alternatives in terms of other kinds of fields. 

Many approaches have been used to reconcile the inherent anisotropies of non-scalar field sources with the isotropy 
of the cosmological background. Ford's original vector inflation [i[ and some more recent vector dark energy models 
[j| rely on fine-tuned initial conditions. In the non minimally coupled vector inflation model, it is supposed that 
there are numerous randomly oriented fields that on the average contribute an isotropic stress energy tensor Q . The 
nonminimal couplings of the model have been claimed to introduce pathologies E l B , see however Q. The model 
l/^ | has also been generalized to nonminimally coupled 2-form and 3- form inflation [lOL lllj . Anisotropies can be avoided 
by setting a "triad" of three vectors pointing in all space- like directions [l2l |. By considering more involved models 
^vq ' with coupled scalar and vector fields, one may also construct dynamical systems where the small anisotropy is kept 
t-H i under control [l^, ■ If the vector field is subdominant, the background isotropy can be retained while the field can 
contribute as a curvaton [l5j] to perturbations in a statistically anisotropic way [l6| ; such models can also be embedded 
in Type II string theory DBI-type inflation [ItJ ■ Time- like vector fields are compatible with the FRW symmetries, 
but a canonical vector is trivial in such a background. Nonminimal couplings and noncanonical kinetic terms have 
• »h j been taken into account in several different context s fl8l - [22j . 

The three-form cosmology proposed in Refs. [TJ l23| can be equivalently described as a noncanonical vector field 
theory. The canonical three-form action can be implemented to successfully model inflation and dark energy (23 - [2fjj 
while respecting the FRW symmetry, as the three-form describes a spatial volume element pointing in the orthogonal 
direction of time. The inflationary models have been found to incorporate large nongaussianities in some parameter 
regions [27j and are claimed to produce more efficient reheating via parametric resonance than usual scalar models 28] . 
Indeed, very distinct features from scalar models arise when any interactions with other fields are taken into account, 
including nonminimal couplings to gravity [ToL ITlj. In particular, a new mechanism for magnetogenesis was proposed 
29] which is based on the simplest U(l) invariant coupling to the electromagnetic field, and this mechanism turns out 
to solve the backreaction problem [30( and to be testable as its predictions for nongaussianity are extremely sensitive 
to the model specifics [31] . Interestingly, such a coupling might be interpreted as an interaction between a hidden 
sector and the visible photons if the three-form was identified as a hidden sector gauge boson in a framework that is 
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very generic in low energy string theory; moreover such coupling can be further constrained by its rich phenomenology 
along the lines of e.g. |33 |. 

In the present paper, we investigate couplings of the three-form to the dark matter particles. Interactions between 
dark energy and dark matter can be justified from high energy physics and phenomenologically they could alleviate 
the coincidence problem. Coupled scalar field quintessence [2J, |3| is very well studied and a myriad of variations have 
been introduced However, in the three- form context qualitatively new effects will appear, due to the fact that the 
coupling not only modifies the conservation laws of the interacting components but will also contribute nontrivially 
to the stress energy tensor, thus modifying the effective source terms for the field equations. Technically this is due to 
the need to contract the three- form indices in order to render the interaction term covariant. This has been neglected 
in all previous investigations of vector-field dependent cosmological couplings [H, H|| and in the three- form models 
of Ref . (35j , since in those studies the interactions were not derived from a covariant Lagrangian but introduced as 
phcnomcnological terms. Such parameterisations have shortcomings: one cannot uniquely deduce the behavior of 
perturbations, and they can even result in unphysical instabilities that have not been found to occur in any proper 
Lagrangian model [36]. What we find here is that in addition to the small-scale instability generic in scalar models 1 
the dark matter will also acquire an effective sound speed due to the novel interaction. 

The contents of the paper are as follows. In the next section |TT] we will present the field equations for the system 
consisting of the three-form field and particles whose mass depends on this field. We model the fluid of particles with 
the point particle action, which is a sufficient description for most cosmological purposes. In section HTT1 we consider 
these field equations in the FKW space-time, and in particular perform a dynamical system analysis of the phase 
that generalizes the previous findings. In section ITVl we derive the perturbation equations for the coupled system and 
analyze it in detail in the Newtonian limit. An evolution equation for the linear structure of matter overdensities 
is a central result of this study. In section [V] we consider some specific example models: first we present a very 
simple case without potential that provides a new kind of scenario in which to address the dark energy problem, 
and secondly we study a model with a nontrivial coupling and a quadratic potential to demonstrate the large-scale 
structure implications in this more generic case. Section [VI] briefly concludes the paper. 



II. FIELD EQUATIONS 

We consider a canonical three-form field ^4 Q ^ 7 with a potential V(A 2 ) coupled to point particles: 

C = -^F 2 - V(A 2 ) £ m a (A 2 )5(x - x(X)), 



(1) 



Here F is the field strength tensor F a pys = ^d[ a Ap^g]- Thus the mass of the particles depends upon the three-form. 
Due to this interaction, an extra term proportional to the coupling will appear in the stress energy tensor, 



T — -F a F aPl + fi 

± iiv g- 1 flap y- 1 v T u 



V'{A 2 ) + 2 Pm (logm(A 2 ))'l A^A vafi - 



±F 2 + V(A 2 ) 



(2) 



where a prime means differentiation with respect to A 2 . We have assumed that the velocities of the particles are 
small, so that their energy density p m corresponds to dust, 



,>„, =Y J ^a{A 2 )8^{x-x{T)) X j S j. 



(3) 



If the extra term is taken to be included within the matter energy stress tensor, one notes it is no longer a perfect 
fluid, however this is a matter of interpretation. The equation of motion for the three form field follows from the 
Lagrangian (TTJ as 



V a F a(3 ~ fS = UU'A^ 5 . 



(4) 



1 This instability, first reported and thoroughly analyzed in ref. [3^ ]. has been dubbed an "adiabatic instability" [38l | (in the sense of 
slowly rolling field), perhaps a bit misleadingly, as it is present in nonadiabatically evolving models as well 13711 and in fact disappears 
in the exact adiabatic limit [39t |40H . This linear instability has been considered also in cases where (part of) dark matter consists of 
neutrinos [4ll I42I ] and was recently shown to exist for generalized nonconformal form of couplings too [43| . It would perhaps be more 
appropriate to call a Laplacian term like we find here an adiabatic instability (in the usual sense of entropy vs. adiabatic perturbations 
in cosmology) since if there were no entropy perturbations, the sound speed squared associated to this term would turn out negative in 
an accelerating cosmology - indeed this is what has proven fatal to some unified models of dark matter [44l |45|| . 
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where we defined 

U' = V'(A 2 )+2p m (logm(A 2 ))' . (5) 

Due to antisymmetry of the three-form, this further implies that 

A a ^U[ a = -U'V a A a ^ . (6) 

Using these equations and the vanishing of the divergence of the total stress energy tensor © then gives us the 
equation of motion for the matter fields: 

V p (p m u"u v ) = ^F a ^ s Qv„F Q/37 * - V a F„p lS ) ~ 2U' (F uaM A a ^ + 3A a ^ a A^) + 2V 1 ' A a ^V a A v ^ . (7) 

Using the antisymmetry of the three-form and its field strength, this further simplifies to 

V M ( Pm u»u v ) = -4 (logm(A 2 ))' A a ^V v A afh . (8) 

This closes the system of equations. It also seems that the coupling vector V M T( mattcr )^ is proportional neither to the 
four velocity of matter u v nor to the four- velocity of the three- form field ~ e-yuB^A^ 1 . This feature is not captured 
by the phenomenological parameterisations considered in the literature, e.g. [36l. kg . 



III. BACKGROUND COSMOLOGY 



The ansatz 



Ajk = a 3 (t)e ijk X(t) 



(9) 



is compatible with the FRW symmetries since it corresponds to the time-like component of the dual vector field. We 
normalize by the cube of the scale factor a(t) so that A 2 = 6X 2 . The potential and the coupling can then be regarded 
as functions of the effective field X. In terms of that field, we define the shorthand notations for convenience: 



/= - (mm) x 

K ' 



\=-QnV) x , 

K • 



where k 2 = 8irG. The Friedmann equations are 



H z = — 
3 



2H + 3H 2 = k 



i (x + zhx) 2 + v + Pn 



\{X + 3Hxf + V-(V X + np m f) X 



x 



The equation of motion for the three-form field is 

X = -3HX - 3HX - V 
and the continuity equation for the matter density is 

p m + 3Hp m = np m fX 



KPmf i 



(10) 

(11) 

(12) 

(13) 
(14) 



As discussed in the introduction, a peculiar feature here is that, in addition to the coupling terms in the conservation 
equations (|13[) and (|14p . there is an extra pressure contribution in the second Friedmann equation (|12p . 



A. Dynamical system analysis 

We define the dimensionless variables 



,Y. ,=^' + 3*). ** = g, «»-^. (15) 
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Here and in the following the prime stands for the derivative wrt e-folding time. In terms of these variables, the 
Fricdmann constraint reads 



1 = y z + z* + uj z . (16) 

This allows us to eliminate one of the variables. We shall consider the phase space spanned by x, y and z. The 
autonomous system of evolution equations can be written as 

x' = -3x + V6y, (17) 

/ (V6-3xy) (y 2 + z 2 - l) + 3y(z 2 (Xx- 1) + l) - V&Xz 2 - 3y 3 , (18) 

-3{fx + 1) (y 2 + z 2 - 1) + 3A.t (z 2 - l) + s/&Xy] . (19) 



y' = l - 

y 2 

z 



Z = 2 

There are four types of fixed points. 

• Type A: Matter domination. Strictly this fixed point exists only when /(0) = 0. It is always unstable since 
the eigenvalues are (|, — ^(3 ± ^81 — 48/ i x(0))). The effective equation of state is of course w e g = for this 
solution. 

• Type B: Three-form saturated de Sitter. This is the cosmological constant which a massless three-form gen- 
erates. Then x = ±^/|, y = 1 and uj = z = 0. The eigenvalues are (—3,-3,0), the vanishing eigenvalue 
corresponding to the perturbations in the direction of z. Plugging perturbations only in this direction and 
expanding to the first nonvanishing order, we get 8z' = —^[3 + \/6(A(-^/§) + f(^J~^))]5z 3 . When the expression 
in the square brackets is positive, this fixed point is an attractor, otherwise a saddle point. The stability thus 
depends upon both the functions / and A. 

• Type C: Potential dominated de Sitter 2 . This fixed point exists when A = 0, i.e. at the extremum of the 
potential at xq. Then y — \J~^Xo and z = y/l — y 2 . The eigenvalues arc (—3, — 1(3± y/9 — 3\,x(xq){2 — 3xq) 2 ). 
In accordance with our intuition, the fixed point is stable if the effective field X has a real mass at the minimum, 
Vxx(x ) > 0. 

• Type D: Scaling solutions. The scaling solution for which u 2 = ||xo//| exists only in the presence of the 
coupling for x\ < |. The universe expands as matter dominated io e fj = and y — \J\x$, z = 0. For this point 
to be stable, we need f{xo)x n < —1. In addition, / cannot be a constant but 

3 (2 - 9x 2 - 4 /( , M3* 2 - 2)) < < 3 / ^18 + 9/(x ) ? + 16\ _ (2Q) 



(2-3x 2 )2 16 V 2-3x 2 

This concludes the quite elegant summary and generalization of all previous results. 

IV. COSMOLOGICAL PERTURBATIONS 

We parametrize the scalar fluctuations of the three- form using the two scalars a and ao, 

Aoji = aajka^, (21) 
A ijk = a 3 e ijk (X + a ). (22) 

The equation of motion, V a F a ^ 5 = 12(V'(A 2 ) + 2p(log m(A 2 ))'A^ 5 , yields the following equation of motion for 

V 2 • • • 

d + 3Ha + (3H + V X x)a a - -^-(a - 2Ha) + (X + 3HX)(3cb - i>) + 3(V XX X - V x )<l> + 2V.x4> = (23) 

- Kfp m (6 m -3<f> + 2ip)- Kf iX p m (oio + 3X(j)), (24) 



a 



2 Note that this coincides with the type B fixed point when z = 0. Also, one could extend the solutions to Xq > f; if negative potentials 
were allowed, corresponding to imaginary z here. 



where S m — Sp m /p rn , and the constraint 



do + 3Ha + (X + 3HX)(30 - V>) + I % - Xr) a = -^^a. 



X a 2 X 



The identity V • (V • F) = gives the additional constraint 



' X - Vxx(a + 3X0) - Vxip = nfp m (S m + 4>) + Kf,xp m {ao + 3X0) ^ 



at \ X J at \ X 

For the matter perturbations we obtain 

L + - - 30 = ^ (/s/(o + 3X0)) , 



and 



^ $ / V 2 

- + — V + fl - = -«/ [X- + — (a + 3X<t>) 



By combining these two we obtain 



5 m + 2H5 rn - 6ff0 - 30 - ^0 = Uh + nfX + -| ( K /(a + 3X0)) 

+ Kf fx(30 - <5 m ) + ^r(a + 3X0) 



For the Einstein tensor we have that 



hG = 3F(0 + #</>) V " 



2 



a 2 

l -5G\ = 0>tf(30 + 0) + (2ij + 3ff 2 )0-i^(0-0), 
2^ = 2?<* 

and the perturbed components of the energy-momentum tensor read as 

- 5T° = Sp = S Pm + (X + 3HX) (a + 3Ha + (X + 3HX){3cf> - ip) - ^-a ) + V,x(a + 3X0), 



5T? = -(Kfp m + Vx)aA + ap m v\ 

5TI = 5p=-{X + 3HX)(a + 3Ha + (X + 3HX)(3(f>-'ip)-^-a) - -V v v A'(n„ - 3A". .) 

nfpm (S m + «o + 3X0) + K,f,xXp m (a + 3X0), 
5T] = 0. 

Einstein's equation 5G\ — k 2 5T^ gives 

= 0. 

Using this result and equations (|3Tj) and (|35|) and eliminating do with ([25]) give, through Einstein's equation 
k 2 6T?, 

V 2 • k 2 ( 6 V 2 

- — (0 + #00 = y ( Pm- - (Kfp m + K^)— a 
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Using flU) and flU} and equating SG^ = k 2 <5T °, yields 



a z 2 



X 



{X + 3HX)a -a - 3X(f> 



! Kj 1 

-p m S m = —nfp m — {X + 3HX)a. 



(40) 



Finally, Eqns. ((MJ) and ((MJ) through 5G\ = k 2 5TI give 



AH<j) + {2H + 3H 

,2 



— (X + 3HX) 



d + 3Ha a T a - 2(X + 3HX) 



2 2 

- 7r^Kxx(ao + 3X(j>) = — —p m {nf + nf tX X)(a + 3X<p) - —nfXp m Sr, 

6 

In the limit of vanishing coupling these equations reduce to those in Ref. [Hj]. 



(41) 



A. The Newtonian limit of perturbations 

Let us then consider Newtonian scales in the late universe. Going now to Fourier space and using (|39|) and (|40|) in 
combination, we obtain in the reference frame comoving with matter, i.e., with 8/a — 0, 



n 2 X 
= - y p m S c m + V x (a + 3Xcf>) -—(V x + Kfp m )a 



(42) 



The comoving density perturbation 8 m ~ S m coincides with that in the Newtonian gauge in the present approximation. 
Constraint (|39|) also tells us that Ha ~ X(j) which means that in the small scale limit, for k/a H, we can neglect 
the last term of (|4"2")l . Keeping this in mind we can use the anti-symmetry condition (f2l)|) to obtain a relation for 
cto + 3X0, 



Q'o + 3X<f) ■■ 



V X X + Kf^xPn 



■0(X<t>, Ha). 



(43) 



Substituting this expression in (|42[) and neglecting the terms proportional to X<p and Ha as they are small when 
compared to the (k/a) 2 4> term for sub-horizon scales, we obtain the modified Poisson equation 



where 



F 



K " (i x V - xF \ A 



V X X + t/,XPm 

Similarly, we can substitute for k/(«o + 3X(f>) « -Fc) m in the equation for the matter density contrast 



2H + nfX 



2F 
1 - F 



1,2 

eff c 2 c 

Om — n PmOm 9 C o£F"« 



(44) 

(45) 
to obtain 

(46) 



We obtained a modified friction term. In addition, the dark matter perturbation feel an effective Newton's constant 
given by the expression 



v eff 



1 -F 



1 



K 2 P 



2 f F+{2H + nfX)F * V ^ jPn 



2 V.XX + Kf,XPn 



Furthermore, the coupling introduces an effective sound speed given by 

F 



-off 



1 - F 



(47) 



(48) 



Recall that the quantity F, which vanishes when the coupling is zero, was defined in (1431) . Due to these modifications 
in the formation of large-scale structure, we expect that the coupling should be quite small in order to produce a 
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viable matter power spectrum. This has been found to be the case for other cosmological models featuring effective 
sound speed terms, such as unified models of dark matter 0, SH| and some modified gravity models [13, IHJ]. If 
the sound speed squared is positive, there will be oscillations at small scales of the matter power spectra, which are 
stringently constrained by observations. We have positive c^ s when 

Thus, there is an instability at small scales unless the slope of the coupling is negative enough. When very tiny, 
such a coupling could be compatible with observed structures and perhaps help to alleviate the possible lack of 
predicted small-scale structure in the ACDM model. The study of specific models and their observational constraints 
from structure formation is outside the scope of the present paper. In the following section we only consider the 
background evolution of the perhaps simplest scenario using a coupling. More general investigation including the 
computation of CMB and LSS will be carried out elsewhere. 



V. EXAMPLE MODELS 



We will look analytically and numerically at two example models, the first specified by a constant coupling / and 
no potential energy, and the other with a coupling and a quadratic potential. 



A. A model with no potential energy 



Here we will present a very simple realization of a three-form accelerated cosmology that may also address some 
of the fine-tuning problems of dark energy. We can accelerate the universe without a potential. In fact, a massless 
three- form is nothing but a cosmological constant 3 . The de Sitter expansion in the asymptotic future is due to the 
cosmological constant term generated by the massless three-form. It automatically adjusts to the current Hubble rate, 
so as long as we get the transition going on around the present epoch, the scale of the effective cosmological constant 
is naturally of the observed magnitude. We consider a scenario where this is achieved by the coupling. 

In the early universe, we set X = X = 0. The coupling will howeve r set t he field moving and eventually it will reach 
the attractor B. Let us consider the simple case m(A 2 ) ~ exp[K/o\/j4 2 /6] so that / = /n is a constant 4 . Setting the 
initial conditions at 1 + z< = 10 10 and fixing / = 1.62 x 10 gives us the realistic evolution depicted in Fig. [T] It 
can be observed that the Universe comes to be dominated by a three-form even though it emerged from a vanishing 
energy density. While the dynamics of the Universe is dominated by a background perfect fluid with equation of 
state parameter 7 — 1 = pb /pb, the equation of motion for X can be rewritten in terms of the dimensionless variable 
y = k{X' + 3X)/V& as 



^'-^ + V^ /e " ; = ' (50) 

where N = ln(a/ao) = — ln(l + prime means differentiation with respect to N, 7 = 4/3 for radiation, 7 = 1 
for dust and we have neglected the contribution of the coupling in the energy density of dust, i.e., we have used 
Pm = Po^m exp(— 3iV). This equation of motion for y has a simple solution 

y = CeW + 2 gm /e 3( 7 -i)„ (51) 
V6(2 — 7) "s 

where C is a constant of integration which is set by the initial conditions. We now take y(Ni) = and match the 
solutions during matter domination and dust domination at N eq = hi(Q r /Q m ). In addition, we consider that px 



3 This fact was used by Hawking in an approach to explain the vanishing of the cosmological constant problem [49| . The canonic term of 
the massless field A is the "gauge fixing" term for the dual vector +A, F 2 (A) ~ (V ■ *A) 2 . Maxwellian electromagnetism extende d by 
adding this term to the fundamental electromagnetic Lagrangian provides an electromagnetic origin for the cosmological constant 50j ; 
if this contribution was excited as a quantum fluctuation during inflation the predicted value could match with the observed one |5l| . 

4 In this case however, the Newtonian limit derived above might not be valid since when defining (145 I t we have assumed that V xx~^f X K P 7^ 
0. 
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5 

log(l + z) 



FIG. 1: Dynamics for a model without potential but with a linear coupling, / = 1.62 x 10 -12 . Evolution of X, y, the energy 
density of the various components of the Universe and the equation of state parameter wx with redshift. In the lower-left panel 
the lines represent the energy density of radiation (dashed line), dust (dot-dashed line) and X (solid line). 



becomes dominant at N a = In^m/flx) 1 ^ 3 , as the equation of state parameter for X is wx ~ —1 at late times. 
Putting all this together we can find a simple estimate for the value of / such that the 3-form field is accelerating the 
Universe today 



3^(1 + Zl )2- 



(52) 



In Fig. [5] we show how this estimate compares with the value obtained by numerically solving the equation of motion 
and demanding fix — 0.7 today. As the coupling parameter /o is very small, the exponential form of the coupling 



10" 




FIG. 2: The solid line represents the value of / required to obtain fix = 0.7 today and the dashed line represents the estimated 
value using Eq. (|52l) . 



is not crucial for our conclusions, but the dynamics would be similar for any coupling that can be approximated by 
m 2 w tuq + JqkA + . . .. 
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B. A model with quadratic potential 

Considering now a quadratic potential for the three form such that V = VqX 2 , the late dynamics of the Universe 
can be fairly different. The presence of the coupling between the three- form and dark energy displaces the field X 
from the minimum towards X w ± y/2/3, driving the Universe to an accelerated expansion. When the energy density 
of the background decays, the field oscillates in the minimum of the potential and the evolution of the Universe once 
again decelerates. The evolutions of X, y, the energy density of the various components of the Universe and the 
equation of state parameter wx with redshift are illustrated in Fig. [3l 




10 5 10 5 

log(l-M) log(l + z) 




10 5 10 5 

log(l + 2) l0g(l+^) 



FIG. 3: Dynamics for a model with a quadratic potential, V = Vo(kX ) 2 with Vo = 4 x 10 _48 GeV 4 , and a linear coupling, 
/ = —1.725 x 10~ 12 . Evolution of X, y, the energy density of the various components of the Universe and the equation of 
state parameter wx with redshift. In the lower- left panel the lines represent the energy density of radiation (dashed line), dust 
(dot-dashed line) and X (solid line). 

When the potential is included, F in Eq. (|45j) is well defined and we can use (|46[) to compute the evolution of the 
matter density contrast. The result is illustrated, for three values of the coupling /, in Fig. 2] as the ratio 5 m (z — 
0)/5m(z = 10 3 ) for scales that at z = 10 3 are k = (10 — 10 3 )aH. The figure shows that there is hardly any difference 
in S m for the various scales, however, below a certain scale, i.e., for large enough k/a, the growth of perturbations 
increases abruptly. To understand this, we point out that F takes negative values and consequently, c 2 s < 0, which 
means that there is an instability in the evolution of the linear density contrast when (k/a) 2 > K 2 s p m S m /2c 2 s . 

In Fig. [5] we show the evolution 8 m (z) with redshift for the same three values of / used in Fig. [4] 

VI. CONCLUSIONS 

We considered coupled three-form cosmologies. In particular, the cosmological fluid was modeled as a gas of point 
particles whose mass depends upon the three-form. We presented the field equations and derived the cosmological 
equations of motion for both the background and perturbations. 

For the background, we performed a phase space analysis generalizing the previous results, finding a new class 
of fixed points and modified stability conditions due to the presence of the coupling. We also presented an example 
model where the massless three-form, initially vanishing, acquires a small, constant energy density as the coupling sets 
it evolving and until the field freezes in the fixed point B. The energy density would match the observed cosmological 
constant provided only that the coupling strength was sufficiently small (around ten-fifteen orders of magnitude below 
the Planck scale depending on when inflation ended). 

For linear perturbations, we analyzed the features at the Newtonian limit and found the evolution of structure to 
be modified by an effective friction term, a time-varying effective gravitational coupling between the particles and 
a sound speed -like effect in the presence of which dark matter is not precisely cold. Since the perturbations are 
so sensitive to coupling, we might see subtle effects even with a very small coupling constant as in the one specific 
model considered here. In particular, the appearance of an effective nonzero sound speed introduces scale-dependence 
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5 m (z = 0)/5 m (z = 10 3 ) 



423 r 



422.5 



422 



421.5 



421 1 



-f = -1.725x10 
-f = -1.725x10 
f = -1.725x10 



-14 



-12 



-10 



10 



10' 



10° 

k/a H 



10° 



FIG. 4: Ratio of the linear density contrast 5 m (z — 0)/8 m {z = 10 3 ) for scales that at z = 10 3 are k = (10 — 10 7 )aH and for 
three values of the coupling /. 



6 m (z)/6 m (z = 10 3 ) 



600 



f = -1.725 x 10 

---f = -1.725 x 10 




2 1.5 1 

log(l + z) 



FIG. 5: Redshift dependence of 8 m for scales that at redshift z = 10 are k — 1.5 x 10 aH for three values of the coupling /. 



to structure formation. The coupling effects will be enhanced at small scales, which in the best case could help to 
alleviate such issues in the small-scale ACDM structure formation as the missing satellite problem. 

Study of large-scale structure constraints in specific models, the workings of possible screening mechanisms and 
disformal couplings j43j remain to be explored in the context of interacting three-form cosmology. 



Acknowledgments 



We thank Anupam Mazumdar and Danielle Wills for useful suggestions. TK is supported by the Norwegian 
Research Council. NJN is supported by a Ciencia 2008 research contract funded by FCT and through the projects 



11 



PEst-OE/FIS/UI2751/2011 and grant CERN/FP/116398/2010. 



[1] E. J. Copeland, M. Sami, and S. Tsujikawa, Int.J.Mod.Phys. D15, 1753 (2006), hep-th/0603057. 
[2] C. Wetterich, Astron.Astrophys. 301, 321 (1995), hep-th/9408025. 
[3] L. Amendola, Phys.Rev. D62, 043511 (2000), astro-ph/9908023. 
[4] L. Ford (1988). 

[5] T. Koivisto and D. F. Mota, Astrophys.J. 679, 1 (2008), 0707.0279. 

[6] A. Golovnev, V. Mukhanov, and V. Vanchurin, JCAP 0806, 009 (2008), 0802.2068. 

[7] B. Himmetoglu, C. R. Contaldi, and M. Peloso, Phys.Rev. D79, 063517 (2009), 0812.1231. 

[8] A. Golovnev, Class.Quant.Grav. 28, 245018 (2011), 1109.4838. 

[9] M. Karciauskas and D. H. Lyth, JCAP 1011, 023 (2010), 1007.1426. 
[10] C. Germani and A. Kehagias, JCAP 0903, 028 (2009), 0902.3667. 
[11] F. S. Koivisto, D. F. Mota, and C. Pitrou, JHEP 0909, 092 (2009), 0903.4158. 
[12] C. Armendariz-Picon, JCAP 0407, 007 (2004), astro-ph/0405267. 

[13] M.-a. Watanabe, S. Kanno, and J. Soda, Phys.Rev.Lett. 102, 191302 (2009), 0902.2833. 
[14] M. Fhorsrud, D. F. Mota, and S. Hervik, JHEP 1210, 066 (2012), 1205.6261. 
[15] K. Dimopoulos, Int.J.Mod.Phys. D21, 1250023 (2012), 1107.2779. 

[16] K. Dimopoulos, M. Karciauskas, D. H. Lyth, and Y. Rodriguez, JCAP 0905, 013 (2009), 0809.1055. 

[17] K. Dimopoulos, D. Wills, and I. Zavala, Nucl.Phys. B868, 120 (2013), 1108.4424. 

[18] C. Eling, F. Jacobson, and D. Mattingly, pp. 163-179 (2004), gr-qc/0410001. 

[19] F. Koivisto and D. F. Mota, JCAP 0808, 021 (2008), 0805.4229. 

[20] F. Zlosnik, P. Ferreira, and G. Starkman, Phys.Rev. D77, 084010 (2008), 0711.0520. 

[21] J. Beltran Jimenez and A. L. Maroto, Phys.Rev. D80, 063512 (2009), 0905.1245. 

[22] J. D. Barrow, M. Fhorsrud, and K. Yamamoto (2012), 1211.5403. 

[23] F. S. Koivisto and N. J. Nunes, Phys.Lett. B685, 105 (2010), 0907.3883. 

[24] T. S. Koivisto and N. J. Nunes, Phys.Rev. D80, 103509 (2009), 0908.0920. 

[25] C. G. Bohmer, N. Chan, and R. Lazkoz, Phys.Lett. B714, 11 (2012), 1111.6247. 

[26] A. De Felice, K. Karwan, and P. Wongjun, Phys.Rev. D85, 123545 (2012), 1202.0896. 

[27] D. J. Mulryne, J. Noller, and N. J. Nunes (2012), 1209.2156. 

[28] A. De Felice, K. Karwan, and P. Wongjun (2012), 1209.5156. 

[29] F. S. Koivisto and F. R. Urban, Phys.Rev. D85, 083508 (2012), 1112.1356. 

[30] V. Demozzi, V. Mukhanov, and H. Rubinstein, JCAP 0908, 025 (2009), 0907.1030. 

[31] F. R. Urban and F. K. Koivisto, JCAP 1209, 025 (2012), 1207.7328. 

[32] S. Abel, M. Goodsell, J. Jaeckel, V. Khoze, and A. Ringwald, JHEP 0807, 124 (2008), 0803.1449. 
[33] H. Wei and R.-G. Cai, Phys.Rev. D73, 083002 (2006), astro-ph/0603052. 
[34] W. Zhao (2008), 0810.5506. 

[35] F. Ngampitipan and P. Wongjun, JCAP 1111, 036 (2011), 1108.0140. 

[36] J. Valiviita, E. Majerotto, and R. Maartens, JCAP 0807, 020 (2008), 0804.0232. 

[37] F. Koivisto, Phys.Rev. D72, 043516 (2005), astro-ph/0504571. 

[38] R. Bean, E. E. Flanagan, and M. Frodden, Phys.Rev. D78, 023009 (2008), 0709.1128. 
[39] P. S. Corasaniti, Phys.Rev. D78, 083538 (2008), 0808.1646. 
[40] P. S. Corasaniti, AIP Conf.Proc. 1241, 797 (2010), 1001.2687. 

[41] N. Afshordi, M. Zaldarriaga, and K. Kohri, Phys.Rev. D72, 065024 (2005), astro-ph/0506663. 
[42] D. Mota, V. Pettorino, G. Robbers, and C. Wetterich, Phys.Lett. B663, 160 (2008), 0802.1515. 
[43] F. S. Koivisto, D. F. Mota, and M. Zumalacarregui (2012), 1205.3167. 

[44] H. Sandvik, M. Fegmark, M. Zaldarriaga, and I. Waga, Phys.Rev. D69, 123524 (2004), astro-ph/0212114. 

[45] F. Koivisto, H. Kurki-Suonio, and F. Ravndal, Phys.Rev. D71, 064027 (2005), astro-ph/0409163. 

[46] L. L. Honorez, B. A. Reid, O. Mena, L. Verde, and R. Jimenez, JCAP 1009, 029 (2010), 1006.0877. 

[47] F. Koivisto, Phys.Rev. D73, 083517 (2006), astro-ph/0602031. 

[48] B. Li, J. D. Barrow, and D. F. Mota, Phys.Rev. D76, 104047 (2007), 0707.2664. 

[49] S. Hawking, Phys.Lett. B134, 403 (1984). 

[50] J. Beltran Jimenez and A. L. Maroto, JCAP 0903, 016 (2009), 0811.0566. 

[51] J. Beltran Jimenez, F. S. Koivisto, A. L. Maroto, and D. F. Mota, JCAP 0910, 029 (2009), 0907.3648. 



